Screw and edge dislocations with time-dependent core width: 
From dynamical core equations to an equation of motion 
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Abstract 

00 . 

' Building on ideas introduced by Eshelby in 1953, and on recent dynamical extensions of the Peierls model for screw 
• ^ . and edge dislocations, an approximate equation of motion (EoM) to govern non-uniform dislocation motion under 
c/2 \ time-vaiying stress is derived, allowing for time variations of the core width. Non-local in time, it accounts for ra- 
T-^ ' diative visco-inertial effects and non-radiative drag. It is completely determined by energy functions computed at 
■ constant velocity. Various limits are examined, including that of vanishing core width. Known results are retrieved 
. as particular- cases. Notably, the EoM reduces to Rosakis's Model I for steady motion [Rosakis, P., 2001. Supersonic 
dislocation kinetics from an augmented Peierls model. Phys. Rev. Lett. 86, 95-98]. The frequency-dependent effective 
response coefficients are obtained within the linearized theory, and the dynamical self-force is studied for abrupt or 

1 . smooth velocity changes accompanied by core variations in the full theory. A quantitative distinction is made between 
' low- and high-acceleration regimes, in relation to occurrence of time-logarithmic behavior. 
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1 Introduction 

rn 

O ' Plastic deformation in crystals arises as dislocations move through the material under an applied stress 
2 (e.g., Hirth and Lothe, 1982). Their individual motion is expected to be determined by an equation of mo- 
tion (EoM), one non-trivial part of which is the self-force (Eshelby, 1953; Ni and Markenscoff, 2008). The 
latter can be seen as the counterpart for defects in crystals of the Newtonian inertial force for particles. 
Dislocation inertia is partly grasped starting from the fact that the field configuration of the material dis- 
5^ \ placement induced by a stationary dislocation is velocity-dependent (Frank, 1949; Eshelby, 1949). As the 
velocity v changes, fields are updated to comply with the new state of motion. This updating proceeds at 
finite wave speed via sound wave emission from the dislocation core, which induces inertia. Updating con- 
tinuously occurs over time in accelerated or decelerated motion, and because a non- supersonic dislocation 
always moves in its own updating wave field, dislocation inertia is history-dependent, i.e., non-local in time 
(Eshelby, 1953; Nabarro, 1967; Hirth and Lothe, 1982, p. 195). Its long-time remanent character, for infinite 
rectilinear dislocations, is due to peculiarities of the associated two-dimensional wave propagation problem 
— the "afterglow" effect discussed by Barton (1989) and Lazar (201 1). 
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The specific question of the self-force and the related question of radiative drag involves dynamical 
fields generated by dislocations. These and related acoustic emission phenomena (Sedgwick, 1968) have 
been considered by a number of authors, notably in connection to strong-motion studies in seismology, 
e.g., Lund (1986) and references therein. To cite but a few analytical works, expressions of displacement, 
strain, or stress fields generated by arbitrarily moving screw or edge dislocations were given for an isotropic 
medium by Kiusalaas and Mura (1964a; 1964b), and recently revisited by Lazar (201 1). Brock (1983, 1986) 
considered non-uniform motion along non-planar paths, in connection with crack extension. Dynamical 
fields produced by a finite dislocation segment have been addressed by Lund (1986). Arias and Lund (1999) 
studied motion in samples of finite width. Beltzer (1982) considered dynamical acoustic emission associated 
to random dislocation motion. The self-force itself has been explored after Eshelby by several authors, 
notably Beltz et al. (1968), and more recently by Markenscoff and co-workers who undertook a systematic 
study of the singularities associated to moving dislocations (Markenscoff, 1980; Ni and Markenscoff, 2008; 
Markenscoff and Huang, 2008). 

Consideration of systems with lattice periodicity leads to involved dispersion relations (e.g., Askar, 
1986), with incidence on dislocation dynamics (Eshelby, 1956). For alternative approaches using gradient 
elasticity see, e.g., Eringen (2002) and references therein, and Lazar (2010). However, as far as inner length 
scales are concerned, lattice-related features can approximately be accounted for in classical continuum 
elasticity via a Peierls-Nabarro (PN) cohesive-zone approach on a prescribed slip plane, which produces 
dislocations of finite width without introducing additional wavemodes. Accordingly, we restrict ourselves 
to planar cores, i.e., Somigliana dislocations (Nabarro, 1967), such as in FCC crystals (Heidenreich and 
Shockley, 1948). 

In this context, this paper addresses the question of obtaining the EoM, with focus on its structure, thus 
trying to answer a question implicitly raised by the phenomenological approach of Pillon et al. (2007). The 
latter suggests that an EoM of simple form, valid for all velocity regimes under minimal assumptions, might 
be within reach. The present work, to be seen as an effort to notably reduce the part of phenomenology of 
the latter work, answers positively. To this aim, we start from dynamical extensions of the well-known static 
PN equation (Peierls, 1940; Nabarro, 1947; Schoeck 2005) for the shape of dislocation cores. Following 
initial steps taken by Eshelby (1953) and, e.g., Stenzel (1969), such dynamical core equations were recently 
obtained for screw and edge dislocations (Pellegrini, 2010; 201 1) by means of the Green's function approach 
(Mura, 1987) in the framework of generalized-function theory. No attempt will be made to solve numerically 
the EoM to be obtained under an applied stress, which is left to future work. We instead try to gain some 
fresh insight into its most important component, the self-force, by artificially prescribing the motion, which 
is usual practice in this kind of studies. The embedding medium is infinite, although finite-size effects may 
notably influence fast dislocation motion (Arias and Lund, 1999; Vandersall and Wirth, 2004). 

The dynamical core equations are reviewed in Section [21 The method used to derive from them the 
EoM is presented in Section [3l Arbitrary time- variation of core width is permitted, a notable difference 
with Eshelby's treatment (1953). In Section|4]the EoM is made fully explicit by means of identities further 
discussed in Appendix A. It involves only known stationary energy functions, which simplifies matters in 
formally bringing inertial effects for edges and screws down to the same level of complexity. The linearized 
theory is considered in Section [51 Section [6l provides further insight in the Volterra limit of vanishing core 
width, in connection with jumps between two arbitrary velocities. Core-width velocity dependence is dis- 
cussed in Section [71 In turn. Section [H addresses velocity changes, either abrupt or smooth in association 
with variations of core width, and analyzes the possible regimes of the self-force. Section [9l specializes the 
discussion to the logarithmic regime. We summarize our findings and conclude in Section [lOl Appendix B 
gathers some useful integrals. Details on numerical calculations are provided in Appendix C. 
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Our conventions are as follows. Shear and longitudinal sound velocities are denoted by cs and cl- A 
screw dislocation is either subsonic (\v\ < Cs) or supersonic (\v\ > cs), whereas an edge dislocation is 
either subsonic (\v\ < Cs), transonic (cs < |f | < Cl) or supersonic (\v\ > Cl). Time and space Fourier 
transforms (FT) are carried out as 

f{x,t)= f^^e^^'^~-y{k,u:), (1) 
J Zn Zn 

f{x, t) being some function of position x and time t. Following usual practice, we distinguish it from its 
FT f{k, cu) only by the symbols of the Fourier momentum k and the angular frequency cu. Integrals with 
implicit bounds run from — oo to +00. 



2 Dynamical core equations 

2.1 Overview 

Let the rectilinear dislocation glide in the plane y = of an isotropic medium of infinite extent. The 
abscissa in the slip direction is x. Let furthermore ?7(x, t) be the relative material displacement (slip) between 
both sides of the slip plane. Given the function aa{x, t) that represents a prescribed resolved shear stress 
applied on the slip plane by external agents, the governing dynamical core equation for ri{x, t) is the stress 
balance equation 

t) + (Td(x, t) + (y^{x, t) = f {r]{x, t)) , (2a) 

where —/'{r]) is the lattice pull-back force that derives from the lattice potential /(?]), usually identified 
to the 7-potential (Christian and Vitek, 1970). This function is 6-periodic, with b the Burgers vector length 
along the slip direction, and such that /'(O) = 0. The function cr^(x, t) represents the dynamical self-stress 
produced by the moving dislocation on the slip plane, and reads (Pellegrini, 2010; 201 1) 

a,{x,t) I dTdx'K{x,t\x',T)^{x',T) - ^^^ix,t), (2b) 

where n is the shear modulus and is a non-local kernel that depends on the dislocation character. Here- 
after, the term proportional to drj/dt, local in space and time, is referred to as the local term. Finally, the 
stress 

a^{x,t)'^-a^^{x,t) (2c) 
Cs at 

has the same form as the local term in (7^,. It represents a phenomenological drag component, with posi- 
tive dimensionless damping coefficient a, which is intended to account for non-radiative damping effects 
(Movchan et al., 1998) and proves necessary for compatibility with a non-zero applied stress in the steady 
state (Rosakis, 2001). We focus on single-dislocation solutions such that ?7(— oo,t) = h + ri{+oo,t) with 
?7(-|-oo,t) such that (Ta(+ 00, t) = /' (?7(+oo, t)), assuming for consistency that cra(+oo, t) = cra(— oo,t). 
We refer to Pellegrini (2010) for the specification of which component of the stress tensor the quantity 
corresponds to, which depends on the dislocation character. 
Expressions for K are as follows. For a screw dislocation. 
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where 9{x) is Heaviside's unit-step function. Instead, for a gliding edge dislocation 
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These causal kernels were derived using a retarded Green function of the material displacement and vanish 
for t < 0. As functions of x, they are regular at x = for t 7^ 0. Because Kl{x, t) is an 0(a;), it is preferable 
in the perspective of integrating over x to emphasize this regularity by writing it in the following form where 
all terms are separately well-defined at x = if t 7^ 0: 
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(5) 



With expressions ^ and ^ of the kernels, Eq. (l2al) has been shown (Pellegrini, 2010) to reduce, in the 
steady-state regime, to Weertman's equationo (1969) in an augmented form with additional drag term. 
The latter extension of the PN model has been introduced by Rosakis (2001) under the name Model I. 



In the above equations, 6'(1 — |x|)/ vT^^a? is merely a particular notation for the composition of 1 



X 



with the locally-integrable function (distribution) x+^^^ equal to x if x > and if x < 0. Consider 
more generally the locally-integrable function x;^" = {x~° if x > 0; if x < 0} where Re a < 1. 

Its derivative is the pseudo-function (^x^"'^ = — aPfx],!""^, where Pf stands for Hadamard's finite part 
(Schwartz, 1966, p. 38). With the shorthand notations 



the kernel K'^ in Eq. (|3]) takes the self-similar form 

e{t) 



K\x,t) 



2cst' 



Ml 



(6) 



(7a) 



where a Heaviside factor 6{t) is introduced to preserve causality in this alternative writing. Moreover, since 
in the sense of distributions 



^ Namely, velocity-dependent stationary generalizations of the PN model. 



4 



Kl{x-x\t-T)— = -^{x-x',t-r) = -^{x-x',t-T), 

by use of (x^"^ = — aPfx^jl""^ with a = 1/2, the kernel K in Eq. (|4a|) is alternatively written as 

K{x, t\x', t) = K^{x — x\t — t) with 

Although we shall proceed otherwise in Section [3] (namely, by using the kernel in its form (|4al)). introducing 
a finite part prescription at this step is quite natural, and makes K"^ well-defined at vs = I- 

Still, due to their self- similarity in x and t, the kernels are ambiguous at (x, t) = (0, 0) where the ratio 
x/t is not defined. The way to go to this limit is a crucial issue in the case of a Volterra dislocation, which 
has no intrinsic length scale. The writings of Eqs. (|7ab) put emphasis on an interpretation of the kernels as 
functions of the "velocities" vs,l defined in and of time. However, this interpretation is possible only 
if t 7^ 0. As a consequence the kernels must be regularized at t = 0. This is conveniently achieved by 
multiplying them by a factor e""/*, where e is a positive constant that sets an "inner" time scale, and by 
taking the limit e — after the time integral over r is done in Eq. (l2bl) . Unless otherwise mentioned, this 
prescription is implicit henceforth; see Pellegrini (201 1) for additional information. 



K'{x,t) 



2c.t2 1 ^,3 



2.2 The local term in 

The complementarity of the regularized non-local kernel K and of the local term in Eq. (l2bl) is illustrated 
in a simple non- stationary case where known results are retrieved. Consider a Volterra screw dislocation 
located at position ^(t), at rest with = for t < and moving arbitrarily at times t > 0. For such a 
dislocation. 



r^{x,t) = he{m-x)- (8) 
Using this expression in (|2bl) with kernel dS]), the stress cr^ becomes such that (Pellegrini, 2010) 
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(t > 0), (9) 



where 



v{x,t,T) 



t-r 



(10) 



The Dirac contribution comes from the local term, and the last term is the integrated contribution of times 
r < prior to motion. 

Now, represent the function x(l — by its Fourier transform — i7rJi(A;), where Ji is the Bessel 

function. With the regularization procedure described above, replacing for convenience e by e/cs, one gets 
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27r 



cs(t-r)2 



(11) 



To make progress, we specialize to the case where the dislocation jumps from rest to a steady state of 
velocity f for t > 0. Then with ^{t) = vt. 



■ k[x — ^{T)]-\-i€ ■ k(x — vt)-\-ie 

t Q cs(t-r) rt g cs(t-r) 

dr — -, = e / dr ■ 



/o"'cs(t-r)2 " Jo cs(t-r)2' 
For x, /c e M, using the change of variables u = l/[cs(t — r)] leads to 



(12) 



■ kx 
g est 



lim / dr — — ■ = lim i- = ie ^=5* p. v. h nS^kx) 

e->0+ JO Cs{t — T)^ e-^0+ KX + 16 KX 



(13) 



where p. v. is the principal value. Since Ji{k)/k is finite at A; = and 5{ax) = 5{x)/\a\ for a G M, a 7^ 0, 
one has 



rt g cs(t-T) 

lim Ji(A;) / dr — - 

e^o+ ^ 'Jo csit-r) 



.Ji{k) / 1\ Ji(A;),, , 

iH^e ^=8* p.v. - + TT— f-^5(x) 



(14) 



so that the limiting value A; = causes no trouble. Applying this result to Eq. (fTTI) . and using the fact that 

Ji{k) I k and J\{k')l\k\ are even and odd, respectively, yields 



lim / dr 



t g cs(t-r) y 
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V y 
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(15) 



The last step consists in making Eq. (1151) explicit by invoking the following integrals (Gradshteyn and 
Ryzhik, 2007): 



; COS{kXj = [1 — X)l. 

IT k 

dk Jijk) 
IT k 



sin(A;x) = x — (x^ — 1)4. sign(x) 



(16a) 
(16b) 



Using the outcome in Eq. ^ and reorganizing terms, the Dirac contribution of the local term eventually 
cancels out thanks to the linear term, x, in Eq. (|16bl) . The self-stress on the slip plane follows as 



27r , , 1 L vt 
—ar,{x, t) = -n + p.v. 

/iO X \ X — vt 



X 



7rsign(t;) (f/cs) — 1 '^^ 6{x — vt). 



(17) 



Apart from an overall minus sign due to a different choice of dislocation sign, the right-hand side (rhs) 
of Eq. (fTTI ) reduces when \v\ < Cs to the stress generated by a subsonic Volterra dislocation. Except for 
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the principal value prescription (Pellegrini, 2011), an equivalent expression is found (up to a factor /i) in 
Eq. (17) of Markenscoff (1980). For > cs we also retrieve Eq. (2) of Callias and Markenscoff (1980) 
at 2 = in their notations. The Dirac term at the dislocation position x = vt, present only for \v\ > cs, 
indicates a coincidence between this position and the tip of the Mach front generated by the dislocation in 
this supersonic regime (Callias and Markenscoff, 1980). 

The cancelation that takes place between the local term and a contribution arising from the integral 
kernels in the above calculation is reexamined in Section |6] from a different point of view, namely, at the 
level of the EoM where the dislocation has a finite core width a. As this introduces an inner time scale in 
the problem, working out explicitly the e-regularization will prove superfluous. 



3 Equation of motion: principle 

An equation of motion for the dislocation position ^ (t) is obtained by eliminating the degrees of freedom 
of the core shape out of the dynamical core equation. To this aim, one multiplies this equation by p(x, t) = 
{dr]/dx){x, t), where ?7(x, t) is its solution, and integrates over x, which results in a work balance equation. 
Because the problem is one-dimensional the latter equation is understood as the EoM for ^(t) (Eshelby, 
1953). Although p is unknown, it is expected that in a first approximation the EoM mainly depends on 
^{t) and on the core width, a{t), rather than of the exact shape of p. As a substitute to the solution t] of 
the core equation, we follow Eshelby and appeal to the usual arctangent ansatz, which solves the steady- 
state equation for a sine pull-back force law (Weertman, 1969), but only for subsonic velocities (Eshelby, 
1956; Weertman, 1967). Accordingly, what follows is not meant to apply to supersonic velocities stricto 
sensu although we shall occasionally refer to the latter to the purpose of enlightening some aspects of the 
equations. 

The Lorentzian p(x, t)/6 can be considered as a delta- sequence to approach the Volterra limit of zero core 
width (Ni and Markenscoff, 2008). Here, however, the value of a is of physical interest. For that reason, we 
relax Eshelby 's original rigid core assumption (Eshelby, 1953) by allowing for a time-dependent core width 
a{t), the latter dependence being left arbitrary for the time being. We thus take 



TT 



2 a{t) 



(18) 



where ^(t) is the position of the dislocation center, and ?7(+oo, t) = rjQ^t). We use henceforth the shorthand 
notation v{t) = ^(t) for the instantaneous velocity. The domain of relevance of this ansatz is not clearly 
established. Numerical comparisons with phase-field calculations (Pillon et al., 2007) have indirectly shown 
that the arctangent approximates well the true time-dependent solution in the subsonic range, for velocities 
less than the Rayleigh velocity cr for edges. No such confirmation is available for cr < \v\ < C]^ and — 
independently of the core shape issue — problems that arise for > cr are evoked in Section |71 in relation 
to the model we shall adopt for a{t). 

Assume that the spatial scale of variation of the driving stress cra(x, t) is much larger than the core size. 
The applied stress averaged over the core, o-^{t), is approximately 

cr,{t) = J dxp{x,t)a,{x,t) ~ a,[^{t),t), (19) 
and 6(Ta(t) is the Peach-Koehler force. On the other hand, the self-force of the dislocation is 



7 



F(t) = I dxp{x,t)aJx,t) = / dT dxdx' p{x,t)K{x,t\x' ,T)p{x' ,t) + -^"^^ (20) 
J TT J 27rcs a{t) 

where ari{x, t) is read in Eq. (l2bl) and where the local contribution, 

--ydxp(x,t)-(x,t) = ^-^ (21) 
has been evaluated using 

^ = - I ^'(t) + k - ^ } P(a^, t) (22) 

and the identities / (\.xp^{x,t) = 6^/[7ra(t)] and / dx[x — ^{t)]p'^{x,t) = 0. The latter three equations are 
consequences of Eq. (fTSi) . By the same token, the phenomenological "viscous" stress in Eq. (l2cl) gives 
rise to a drag force 

Fo t = 2«f ^. (23) 

Due to relativistic effects a depends in general on the velocity, so that this drag force is nonlinear in v (see 
Section IT]). 

Finally, because / is 6-periodic and ?7(— oo, t) = b + ri{+oo, t), 

J dx p{x, t)f'{r]{x, t)) = fivi+oo, t)) - f{r]{-oo, t)) = 0. (24) 

Thus, the core equation projected on p becomes the force balance equation 

F{t) + Fu{t) = ba,{t). (25) 

Since F(t) is a functional of ^ over its history, this equation constitutes the desired EoM for ^. The self- 
force is tantamount to a visco-inertial drag force (Nabarro, 1951; Eshelby, 1953). The following focuses on 
obtaining its expression in a physically appealing form, and on examining some consequences. 

For practical purposes, it is useful to transform the double integral over x and x' in Eq. (l20l) in the 
following way. Let p{k, t) be the spatial FT of t). Setting Ax = x — x' and At = t - r > 0, one has 
for a translation-invariant kernel Ki of the or types 

j dxdx p{x,t) Ki{x-x',t-T)p{x\T) = J dAxKi{Ax,At) J ^ p{k,t)p{-k,T)e'''^\ (26) 
Introduce now the average complex velocity between instants t and r defined as 

^(t,r)^'^^^^^ + i^, (27) 
^ ' ^ At At ' 

where a(t, r) = [a{t) + a(r)]/2. With the ansatz (fTSi) . one has 

p(A;,t) = -6e"'«W^"5'^(*)l'=l. (28) 
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Since Imv > 0, one finds 



— p{k, t)p{-k, T)e"'^' = - Re dk e-(^-^*^)'= = Re ^ ^ ^ 

271 



TT Jo 



At 



whence 



J dx dx' p{x, t) Ki{x — x',t — t)p(x', r) = b^Re J d 

In a similar way, since d/dx' = —d/dAx, one has 

J dx dx' p{x,t)K^{x - x' ,t - t)^{x' ,t) = -^^^ J ^ 



Ax Ki{Ax,At) 



ITT 



Ax K^{Ax,At) 



At 



ITT 



(29) 



(30) 



(31) 



Roughly, these integrals over Ax/ At perform the task of replacing within Ki (resp., in dK2/dx) the quan- 
tity Ax by v(t, r) At. We now are in position to compute F(t) for arbitrarily moving dislocations. 



4 Self-force for screw and edge dislocations 

4.1 Screw dislocation 

Invoking Eq. (l30l) with kernel ([3]) appropriate to a screw dislocation, and appealing to (IB.lal) to carry 
out the integration over Ax/ At by changing variables with u = Ax / (csAt) , the generic expression (l20l) of 
F{t) takes the form 

^ ^ 27rcs i-oo At2 /i_^2/^2 27rcs a{t) 



where v stands for v{t, r). Introducing the characteristic energy per unit length of dislocation line as (e.g., 
Hirth et al., 1998) 

wo = /i&V(4vr), (33) 
it is realized at this point that the function 

Wo (v/cs) 

viy) = , (34) 



is nothing but the quasimomentump = dL/dv, where L{v) is the screw Lagrangian function (see Appendix 
A). Thus, 
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Introducing moreover the associated stationary mass function m{v) as 



def dp 

m = —, (36) 

dv 

one partial integration applied to the integral yields 

Indeed, provided that the dislocation was stationary in the remote past, with 'initial' velocity Vi, the bound- 
ary term Re\p{v) / At] in the partial integration vanishes at r = — oo where v{t, t) = Vi + iO"^. Its vanishing 
at r = t~ as well stems from the following considerations. Eq. (ITTI) entails the expansion 



v{t, t) = i-^ + lit) - U{t) + 0(r - t) (38) 
t — T 2 

from which one deduces, with the asymptotic expansion (IA.3I) . that 

p{v) Wo i 



At cs At 



0(At). (39) 



The real part of this quantity goes to zero as r — t~, as announced. We thus reach our main result: with the 
arctan ansatz and given any a priori time dependence of the core width, the self-force can be written in the 
mass form 

FW = 2Re/' ^rnm,r))f(t,r) + 2^^% (40) 

j-oo t — r dr Cs a[t) 

Thus, the dynamic self-force is expressed exclusively in terms of a mass function already known from sta- 
tionary calculations. Formerly, Beltz et al. (1968) addressed the dynamical problem in such terms, stating 
that "the radiation output of an oscillating dislocation is entirely determined by the dislocation mass factor 
m". The present equation shows that this statement can be generalized. However, Eq. (l40l) . while resem- 
bling their Eq. (37), differs from it in two major ways: the function v(t, r) enters the integral, rather than the 
instantaneous velocity f (r), and the upper limit of the integral in Eq. (|40|) does not need any phenomeno- 
logical cut-off, in close connection to the presence of the rightmost additive term (see Section [6]). 



4.2 Edge dislocation 



The edge case is addressed in the same way. The appropriate kernel is given by Eq. dH), with /i| as in 
(|4c1 ). and Kl as in ([5]). To compute the integral over x and x' in the generic expression of F(t), Eq. (|20l ). use 
is made of Eq. (|30l ) for the contribution generated by Kl, and of Eq. (|3T| ) for that generated by A'g. Consider 
first the contribution of Kl, that comprises four parts. Each of these corresponds to one the four main terms 
in Kl{x,t), see Eq. which are integrated by appealing, respectively, to integrals (IB.lcl) (used twice, 
once with cs and once with cl), (IB.ldl) and (IB.lal) after elementary changes of variables. The contribution 
of Ji| is integrated using (IB.lbl) . Gathering all contributions, immediate simplifications yield F(t) in the 
form (1351 ), but now with 
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p{v) i^^fCs^^ 



v\ 2 - {v/csf 



CsJ [l-iv/csrf 

(41) 



Again, p = dL/dv with the function L appropriate to edges (Appendix A). The steps leading to Eq. MOl) 
can be reproduced, using the asymptotic expansion (IA.3|) to justify the vanishing of the boundary terms in 
the partial integration. In this way, the self-force for edges is cast in mass form as well. 

4.3 Comments 

First, it is emphasized that the mass functions m{v) that appear in the problem are identical to those 
given by Hirth, Zbib and Lothe (1998), but with logarithmic factors removed, as discussed in Appendix lAl 
Also, the kernels K^'^(x, t) in Eqs. (fTab ) stand as close relatives to the quasimomentum p{v) in Eqs. (|34l) . 
(|4T]) . Indeed, 

K(x,t) = ^rVfTV (42) 



2wq \ t 

In this writing, p{x/t) must be interpreted as a distribution with, in the edge case, a finite part prescription, 
in view of Eq. (iTbl) . 

The structure of F(t) in mass form can be understood as follows. In Eq. (l40l) . the functions m{v) involve 
terms of the type (1 — v'^/c'^Y^'^, in which the velocity v is compared with a wave velocity c to determine 
its sub- or supersonic character with respect to c. Heuristically ignoring the presence in v of the imaginary 
contribution of a, the replacement in m{v) of v by 

(43) 

accounts for the fact that whether the dislocation at space-time location {^{t),t) goes faster than a wave 
emitted at (r, ^(r)) for r < t, or not, depends on the average velocity of the dislocation between these 
instants rather than on the instantaneous velocity v(t). Of course, v{t, r) tends to v{t) in the limit of vanish- 
ing time intervals. Moreover (again heuristically), the time-integral in (|40|) is tantamount to a sum of "mass 
times acceleration" terms, where the role of the acceleration is played by 2dv/dT. The factor 2 in the latter 
is necessary to recover the instantaneous acceleration at small times (Pillon et al., 2007). Indeed, as r — t~ 
and again ignoring a, one has 

2*(, r) . 2 ^"' = m + ^mr - <> + o ((. - tf) . m 

The logarithmic integration over past times can be justified from dimensional considerations. Finally, the 
local term proportional to v{t)/a{t) in the self-force regularizes a divergence in the integral as a — in the 
relaxation regime at large times, as will be made clear in Section [6l 

4.4 Link with a result ofEshelby (screw dislocation) 

As recalled in the Introduction, besides proposing his dynamical core equation for screws, Eshelby 
(1953) proposed in the same paper an independent derivation of the self-force of a screw, from an elec- 
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tromagnetic analogy between a moving screw dislocation and a moving current line (his Eq. (26) in that 
reference). He used it as a starting point for further analysis. This equation relies on the arctangent ansatz 
with a rigid core a{t) = a = const. Under this constraint, we show that Eq. (|40l) matches Eshelby's. In our 
notations, and with = ^(t) — ^(r) and At = t — r, Eshelby's self-force F{t) read^ 



m 



2 Re 



dr 



dv 
d7 



(r) +4- virf 



d 



(45) 



where we have carried out the integral over Fourier wavevectors in his Eq. (26). Using integration by parts, 
the dislocation being at rest in the remote past. 



F(t) = 2^Rer dr A 
Cs J-oo \ dr 

d 



Cs 
+ 



v{t) 



viry 



clAt^ - (A^ + ia)2 
d 1 



Me "'^"^dr 



Cs a 



2woCs Re / dr 



4At^ - (A^ + ia)2 

A^ + ia - v{T)At 
[4At^ - (Ae + ia)2]3/2' 



(46) 



where v{t)/a stems from the boundary term at r = t, and where the remaining term under the integral 
has been simplified by evaluating the derivatives. Specializing expression (|27l) of U to a rigid core, namely, 
writing 



v{t, r) = {Ail At) + i(a/At) 



(47) 



and observing that in this case (df/dr)(t,r) = \v{t,T) — t>(r)]/At, Eqs. (140)) and (l46l) are seen to be 
identical. 



5 Linearized equation of motion 

For small velocities, namely, when both v (t) and A^ / At = Re v{t, r) are small compared with a{t, r) / At 
and with the sound velocities, the self-force is linearized by expanding it to first order in these quantities, 
considered as of the same order. The EoM is then solvable in terms of Fourier or Laplace transforms (Es- 
helby, 1953; Al'shitz et al., 1971). Using definition (ITTI) of v{t, r), we assume that the dislocation width 
depends on time only through the velocity: a{t) = a{v{t)). This assumption is discussed in Section |7] where 
the function a{v) is determined. 

Here and in the rest of this section a stands for 5(0), the core width at rest. To linear order in A^/At and 
V, one has 



Re 



m{v) dv 1 
At dr ~ At 



a \ d A^ a , / a \ A^ 
AtJ dr At At2 V At J At 



d 

dr 



m 1 



a\ A^ 



At J At2 



— mi 



a \ Ae 



At; At3^ 



(48) 



Expressions of quite a similar structure were obtained by Lazar (2010), in the context of gradient elasticity. 
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where, since m{v) essentially depends on v"^, m{ia/ At) is purely real. This equation makes clear that the 
time-dependence of a{t) plays no part to linear order. Because the boundary term vanishes upon integrating 
( |48] ) over r the linearized self-force reads 



F'^-{t)=2—-^ -2 drm i— 



a Cs J-oo \ AtJ At^ 

2^ew-2^ew + 2r ^"6^1 ?(^). (49) 



acs J-oo At'^ \ At 

In going to the third line, a change of variable u = ia/ At has been used, together with the identity m{v) = 
W'{v)/v between m{v) and the total line energy W{v). This relationship is discussed in Appendix A where 
W{Q) is given and where it is shown that W{ioo) = 0. The result is of the form 



F'-(t) = -Com + cm + dr C^it - r)e(r) 

J —00 



(50) 



where Co, Ci are constants, and where C2 is a causal response kernel with C2{t) = if t < 0. Fourier- 
transforming with respect to time yields the damped linear oscillator form (Nabarro, 1951) 

F""(a;) = [-io^m{uj) - iioa{uj)]m). (51) 

where rh(uj) = [Co — ReC2{uj)]/uj'^ is the frequency-dependent mass and a{uj) = Ci + liaC2{uj) /uj is the 
frequency-dependent damping coefficient, two real quantities for to real. The tilde distinguishes them from 
the previously introduced mass m{v) and phenomenological damping constant a. 

For the screw, m(a;) anda(u;) were exploited by Al'shitz et al. (1971) on the basis of Eshelby's linearized 
self-force, Eq. (28) in (Eshelby, 1953). They were revisited by Pillon et al. (2007) who spotted an incorrect 
1/2 factor in the linearized expressionHknd started from Eshelby's Eq. (26) instead. Here, this amounts 
to proceeding from the equivalent Eqs. (l40l) or (|46|) . Introduce characteristic times ts,L = a/cs,L' ^iid the 
reference mass 

mo = Wo/cl= ;u6V(4vrcs). (52) 
Equation (3) in Pillon et al. (2007), namely. 



is nothing but the particularization of Eq. (l49l) to screw dislocations. Let /„,, Kn and L„ denote the modified 
Bessel and Struve function (e.g., Abramowitz and Stegun, 1972), and introduce for convenience auxiliary 
functions IL,, as 



This explains the discrepancy between Eq. (I58ab below, and Eshelby's (1953) Eq. (30). 
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(z/ = 0,l,2). 

The following expressions were obtained by Pillon et al. (2007) from Eq. (l53l) : 



~ — log 



1 - \uj\tsKii\uj\ts) 



,2e2-T 

-2 



0(u;' logical), 



2(tsw) + exponentially decaying terms, 



(54) 



(55a) 

(55b) 
(55c) 



mo 



2 1 



TT 



/i(|a;|ts)-L_i(|a;|ts) \ = 2ILi{\uj\h) 



IT 



2 + 0(cj^ 



(55d) 

(55e) 
(55f) 



(7 is Euler's constant). The identity lii{z) = L_i{z) — Ij-n has been used in Eq. (|55d| ) to further reduce 
5(0;). Equations (I55al) and (I55dl) are readily recovered upon Fourier transforming Eq. (|49l ) with Eq. (lA.lcl) 
by means of integral (IB.2bl) . 

To our knowledge, the frequency-dependent response coefficients of the edge are not available (see 
however Kiusalaas and Mura, 1 964b). The Fourier transform of Eq. (|49|) combined with (IA.2cl) is carried out 
with the help of integrals (IB.2al) - (IB.2dl) . After some straightforward simplifications involving the identities 
Kq{x) = -2Ki{x)/x + K2{x), loix) = 2h{x)/x + hix) and Lo(x) = 2'Li{x)/x + L2(a;) + 2x/(37r) 
(Abramowitz and Stegun, 1972) one arrives at: 



m[u) 



mo 



coHH 



1 + 2\uj\tsK,{\u;\ts) + 6K2{\u;\ts] 



/2 



log 



1 + 2\u\ti^K,{\uj\ti^) + 6K2i\co\t^) 



2Ko{\co\ts) 



0(0;^ log 



t 



4 l-|j(ts.)' 



exponentially decaying terms, 



(56a) 
(56b) 
(56c) 



aiu) 



mo 



6/L2(|c^|ts) -2|c^|tsiLi(|a;|ts) 



6iL2(|a;|tL) - 2|a;|tL iLi(|a;|tL) + 2|a;|tsiLo(|a;|ts) 



f^s 

|(l + |)|.fe + 0(.' 
2 + 0(u;-2). 



(56d) 

(56e) 
(56f) 



14 



a 




-0.8- 



b 




Fig. 1. Top: normalized frequency-dependent effective mass fh (solid, black) and damping coefficient a (solid, grey) 
for (a) screw, and (b) edge dislocation for ti^/ts = 0.7. The low-frequency and asymptotic behaviors expressed by Eqs. 
(I55bl ). (I55cl ). (I55el ). and (I55fl ) for the screw, and corresponding equations for the edge, are dashed. Bottom: relative 
difference between screw and edge for (c) m{u)) and (d) a{uj), plotted vs. frequency for ti^/ts = 0.0, 0.3, 0.48, 0.58 
and 0.7 (line thickness increases with tL/ts)- The range of interest for most metals (0.48 < t^/ts < 0.58) is shaded. 

Equation (I56bl) is rewritten as 



-m(0)log(|a;|r 



t* 



9, V 



1 



1 

l+<# 



1_ 



LI 



(57) 



where cf) = {ti^/ts)"^ and m{0) = (1 + (p)mQ is the rest mass factor of the edge. 

By employing Laplace transforms (LT) of variable s instead of Fourier transforms, Eq. (ISTI) becomes 
= s^m(s)^(s). In view of the relationship between FTs and LTs the function m is such that 
m(— icj) = rn(uj) + ia{uj)/uj. The FTs are then deduced from 



m(s) 
mo 



[stsBTiists) 



(screw), 



(58a) 



''S 



6Hr2(sts)-2stsHri(sts)-l 
6Hr2(stL)-2stLHri(stL)-l 



2Hro(sts) (edge), 



(58b) 



where the auxiliary function 'Hy^{z) = | [Hi,(2;) — V^l^)] is defined in terms of the Struve function Ji^, 
and Bessel function of the second kind (Abramowitz and Stegun, 1972). The function m(s) is free of 
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singularities in the complex plane, except for a cut on the negative real axis. It will be used in Section |9l 

The functions m and a are drawn in Figs.[TJa), (b) for ti^/ts — 0.70, i.e., Cl — 1-4 Cs, a situation where 
the screw and the edge behave in quite a similar way. The damping coefficient a{uj) is always positive. 
So is m'"'=*(w). The ratio tt/^s = Cs/cl = [2(1 - z/)/(l - 2i/)]-i/2 decreases from 1/^2 to when the 
Poisson ratio u increases from to 1/2, so that cs < cl/v^. In this range, rff^^^{uj) > 0, as illustrated 
by the asymptotic value in Eq. (I56cl) . which changes sign at cs = cl. For metals varies between 0.032 
(beryllium, ti^/ts — 0.69) and 0.45 (thallium ^l/^s — 0.30), while it stands in the range 0.25-0.35 for most 
of them (Lide, 2005). This corresponds to ^l/^s ~ 0.48-0.58. Comparisons between screw and edge for 
m and a are displayed in Figs.[2c) and[TJd) for < ^l/^s ^ 0.7. Positive values indicate that the screw 
dominates over the edge. While various situations are met for m{uj), depending on the frequency range and 
on the Poisson ratio, edge dislocations experience in general more radiative damping than screws at low 
frequencies. The situation is reversed at high frequencies. These observations remain somewhat formal, in 
view of our using planar cores and isotropic elasticity. 



6 The local term as a counter-term 

The purpose of this section is to show that the local term in Eq. i40\i . inversely proportional to the 
dislocation width a, acts as a counter-tern\^ (CT), which compensates for one specific contribution to the 
self-force integral that explodes at finite times as in the Volterra limit a — )■ 0. This property relates to 
the observations made in Section [Z2l while computing the stress cr^ of a Volterra screw. 

Although we derived it for an isotropic medium, let us assume that the expression of the self-force in 
its mass form is of general validity, and consider as known the total energy function W{v) (see Appendix 
A). The game to be played is to pretend that that the local contribution is unknown. Writing it as CT, the 
self-force reads 



The unknown CT stems from requiring that an expansion of F(t) in powers of a does not lead to infinities 
as the stationary limit is approached in the subsonic regime. 

Indeed, let the dislocation undergo at t = an instantaneous velocity jump between arbitrary initial and 
final velocities Vi and Vf, while it passes at the origin, so that ^(0) = 0. We use the symbol |/] = /f — fi 
to denote the jump of any quantity / between some initial and some final state. We take v{t) = Vi 9{—t) + 
Vf6(t), and suppose first for simplicity that a is constant. This assumption is relaxed in the next sections. 
Then, with \v\ = Vf — and for positive times. 




dr 



dv 



(59) 




if 



if 



< r < t. 



r < 0, 



(60) 



We carry out the change of variables u 



v{t, r), such that du = (dt'/dr)dr. Thus, 



^ The denomination is boiTowed from quantum electrodynamics. 



16 



t-T 



if r G (-00, 0], i.e., for u on a path [vi, Vf + i(a/t)]; 
if r G [0, t], i.e., for m on a path [vf + i(a/t), +ioo). 



(61) 



Using identity (IA.8I) . this leads to the following expression of the self-force: 

"f+if du{u - Vi) dW 



F - CT = 2 Re 



2 Re 



(|f ] t + ia)u du 

W{vf + i^)-W{Vi) 



du — [u) 

lau du 



M t + la 



\v\ t + ia 



'^f+if du dW 



+ 



W{i00) -Wivf + i^j rioo dwdW^ 



la 



ia Jv+i^ u du 



u du 



u 



u) 



(62) 



The remaining integrals are readily evaluated in terms of the kinetic energy function of the field, W]^, 
thanks to Eq. (lA.lOl) . and further simplified in terms of the "stationary Lagrangian" L{v) (|A.4I) . and of the 
momentum relationship p{v) = 2W\^{v)/v, see Eq. (|A.6I) . One ends up with 



ACT 

a M + if 

ACT CT - ^ Im [v,p{u) - W{u)l^,^ . 

(Ji 



(63a) 
(63b) 



The self- force is now straightforwardly expanded to leading orders in powers of a/t, and simplified using 
dL/dv = p. This expansion, which probes the relaxation regime after the jump, is relevant also in the 
Volterra limit a — 0. It reads 



2 2 
F = ACT + - Im L{vi) + - Re 
a t 



p{Vi) 



m 



0(a). 



(64) 



For it to provide a finite leading-order term for \vi^f\ < cs when all energies are real, it is necessary that 
ACT = 0, from which the counter-term must be 



CT = - Im (t)p(ioo) - W{ioo)] , 



(65) 



where v(t) stands for Vf since t > 0. It is shown in Appendix A that, for both screw and edge, W{ioo) = 
andp(ioo) = 2[W]i{u) / u]u=ioo = iuJo/cs- Thus, the local term of EoM is retrieved exactly, namely, CT = 
2wov{t) I (acs). This calculation suggests that its general expression is 



CT = 2Im[p(ioo)] 



vitl 
aity 



(66) 



Eventually, the asymptotic self-force is given by Eq. (l64l) with ACT = 0, which is implied in the rest of the 
paper when referring to (l64l) . The way it is obtained indicates that fi,f should be interpreted as fi,f + iO+, 



17 



where the infinitesimal positive imaginary part is required to get correct determinations of the energies 
above Cs. 

For velocities \vf\ > cs, the term (2/a) ImL(t>f + iO+) in Eq. (|64l ) is non-zero and represents the drag 
force due to concentrated radiation at Mach fronts (Weertman, 1969), to be counter-balanced in steady 
motion by the applied force (e.g., Rosakis, 2001). Its blowing-up as a —t- - due to a constant amount of 
energy on an infinitely thin front - together with that of the phenomenological drag force (l23l) . suggests 
that any kind of instantaneous drag stress not proportional to the dislocation width or to higher powers 
thereof leads in the limit to an infinite self-force contribution, consistently with the presence of the Dirac 
term in Eq. (fTTl ). 

Moreover, bearing in mind that p = dL/dv, the quantity |L] / |y] represents the quasimomentum of the 
dislocation a? jump time, which stands as the relevant 'initial' impulsion prior to relaxation. Consequently, 
p{vf) — [L] / |f ] in Eq. (|64|) is the remainder of the impulsion transfer. We emphasize that Lagrangians 
with imaginary parts are only rarely encountered in classical mechanics, one noteworthy instance being in 
Dekker's (1975) approach to dissipative systems by a complex Lagrangian. 

Equation ( |64l ) embodies former results by Clifon and Markenscoff (1981) for subsonic dislocations 
jumping instantaneously from rest {vi = 0) to a velocity v = Vf. Indeed, since Wk(0) = 0, one has 
L(0) = -W{0) and L{v) - L(0) = vp{v) - [W{v) - W{0)], so thatEq. ^ reduces to 

F{t) = 2K{v)/{vt), (67) 

where K{v) = W{v) — W{0) is the kinetic energy (lA.l II ) of the stationary dislocation. By using expressions 
(lA.lcl) and (IA.2cl) of W{v) in Eq. (l67l) and reorganizing terms, one sees that Eq. (l67l) matches Eqs. (40) and 
(42) (screws), and (47) and (49) (edges) in the above reference. Wu (2002) provides a generalization of Eq. 
(|67]) to anisotropic media. 

Finally, it can be instructive to retrieve this Eq. (|67] ) directly from Eq. (fTTl ) (screw case). Let x = vt + ein 
ari{x, t), and carry out a power expansion of in e. The first-order term is proportional to 1/e and represents 
a divergent contribution to be ignored, thanks to the principal value prescription. The next term is a finite 
constant. Multiplying it by p{x, t) = —b5{x — vt) and integrating over x yields Eq. (|67l) back. 

7 Stationary approximation for the time-dependent core width 

Before going on with dynamical regime changes, we must discuss the time-dependence of the core 
width. To obtain the EoM from the core equation, a "projection" of the latter on p{x, t) was used in Section 
|3l conjointly with an arctangent ansatz of variable width a{t). Although this would be required to complete 
the theory, we shall not attempt to derive an evolution equation for a{t) by like means. Instead, we continue 
to assume an instantaneous dependence a(t) = a{y(t)), written as a{v) for brevity. Some limitations of this 
assumption are however pointed out. 

This assumption is the "minimal" one leading to a meaningful stationary regime. On the one hand "rel- 
ativistic" core distortion evidently imposes some well-defined dependence on velocity; on the other hand, 
the issue of the stationary EoM has been examined by Rosakis (2001), who pointed out that the stationary 
PN equation being conservative, stress-driven subsonic stationary motion is impossible without adding in 
some phenomenological friction term. This resulted in his Model /, already referred to in Section [2?T1 which 
fully determines a{y), friction included. 

The friction stress (l2cl) gives rise to the drag force Fd (l23l) in the EoM (l25l) . The latter is consistent 
with Model I at stationarity: the stationary self-force results from letting t — oo in Eq. (|63al ) or more 
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simply in Eq. (l64l) . which is also relevant to large-time behavior. Accounting for Fd yields the stationary 
stress/velocity relationship 



4a^— + ^lmL(v + iO+) = ba^. (68) 
a{v) Cs a{vj 

To make the connection with Model I conspicuous, introduce the theoretical shear strength (Jth = /ife/ {271(1) 
(Hirth and Lothe, 1982) where d is the interatomic plane separation, so that 2wq = dbath. Introduce next 
functions C{v) (complex), and A{v) and Ba{v) (real) defined as 

C{v) = -A{v) + iB^iv) = -^L{v + iO+) + ia-, (69) 

2wo Cs 

in terms of which Eq. (|68l) is rewritten as 2[d/a{v)]Ba{v) = ojo^a- Introduce moreover the modulus 
D{v) = \C{v)\ = [A\v) + 52(w)]V2. Finally, define a function Biv) = Bo{v) for a = 0, such that 
Ba{v) = B{v) + a;(i;/cs). For c > 0, and with the principal determination of the square root, one has for 
real t;, as e — )■ 0+, 

- {v/c + iey = ^J\l-v^/c^\ [9{c - \v\) - i sign{v)9{\v\ - c)] . (70) 

From this, and by specializing L{v) to screw and edge dislocations using Eqs. (IA.5h b). one verifies that the 
functions A(v), B(v), Ba{v) and -D(f ) in Eq. (|69l ) are exactly those in (Rosakis, 2001), with A(v) and B(v) 
acting as coefficients of the nonlocal and local terms in the augmented Weertman equation of Model I. The 
function A(v) is nonzero only for \v\ < cs (screws) or < cl (edges), whereas B{v) is non-zero only for 
|w| > Cs (screws and edges). 

In this framework, the function a(v) is readily obtained. As mentioned in Section flAl the dynamical 
core equation (l2al) reduces to Model I at stationarity for a sine puUback force fit]) = (Xth sin(27r?7/6). Then, 
a necessary condition for the arctan ansatz (fTSi) to be a solution of the latter is that a{v) = 2dD{v), where 
the above-defined D function is proportional to the energy dissipation rate (Rosakis, 2001). This "external" 
argument is our present substitute for the missing evolution equation for a{t). With this a{v), Eq. (l68l) 
reduces to Rosakis's kinetic equation: 

B^iv)/D{v) = aja,^ (< 1). (71) 

The inequality stems from the relationshipbetween A, Ba and D and indicates the breakdown of the 
single-dislocation solution when > athliJ From the above, the core width and the kinetic relation are 
expressed in terms of the Lagrangian and parameter a in a succinct (but equivalent) formulation alternative 
to Rosakis's, 

a{v) = 2d\L{v + iO^)/{2wo) +ia{v/cs)\ , (72a) 
da = a{v) = ath sin Arg L{v + iO'^)/{2wo) + ia{v/cs) , (72b) 

with, in particular, a(0) equal to d (screw) or 2d[l — (cs/cl)^] = d/{l — u) (edge). 



^ With drag term added, Weertman's equation admits then a "staircase-like" solution, which represents a train of 
kinematically-nucleated supersonic dislocations. Briefly evoked by Eshelby (1956), it has been explicitly given by 
Movchan, BuUough and Willis (1998) in a different context. 
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a = 0.2 




Fig. 2. Core width a{v) vs. dislocation velocity v in Rosakis's Model 1, Eq. ( |72a| ). for a = and a = 0.2. Black: edge 
dislocation; grey: screw dislocation. Shaded zone: a-dependent instability domain v G (cr, c*). A value cl = 2cs is 
used. 



Graphs for cr(t>) are given by Rosakis (2001). The supersonic regime necessarily takes place at saturated 
stress value a{v) = athi^ Fig. [2] represents the core width a{v), Eq. (|72a| ). for drag coefficients a = and 
a = 0.2 (see also Fig. 1 in Pillon et al., 2007). In the absence of drag, a{v) vanishes either at w = cs (screw) 
or at the Rayleigh wave velocity cr ~ 0.93cs (edge), cr being the solution v of A{v) = (Eshelby, 1949). 
With drag added, a(v) is everywhere strictly positive. 

Fig. [2l illustrates the fact that the edge core width blows up at = cs, which dramatically reduces the 
radiative drag, and points towards some pathology of the arctangent ansatz- This blowing-up probably con- 
stitutes the signature of a rich dynamical process involving core dissociation into several partials, followed 
by recombination as Cs is overcome, as reported in some molecular dynamics simulations (Li and Shi, 2002; 
Olmsted et al., 2005). Obviously, a single-dislocation ansatz cannot be much accurate in trying to capture 
such an event, but the mere existence of some analytical signature is worth mentioning. 

However, within Model I for an edge, an instability range Cr < |f | < c* exists, identified by a negative 
derivative da/dv < 0, where cs < c* < cl is the solution of da{v = c^)/dv = (Rosakis, 2001). The 
upper bound c* depends on a, and decreases from the "radiation-free" velocity c* = a/2cs (Eshelby, 1949; 
Weertman, 1969; Gumbsch and Gao, 1999; Gao et al., 1999) for a = down to Cs as a increases. No steady- 
state is possible for such velocities, which probably makes inadequate the approximation a{t) = a{v{t)) 
there. Since this range overlaps the interval [cr, cl], and for lack of a dynamical governing equation for a{t), 
we shall only consider velocities less than cr in numerical applications. This issue concerns a{t), but not 
the foregoing calculations. 



° This unrealistic feature of Model I is alleviated in Rosakis's Model II with gradient (2001), which admits larger 
sti'esses. No such coiTcction is attempted here. 



20 



8 Velocity changes with varying core width 



8. 1 Simultaneous jump in velocity and core width 

The calculation of Section [6] extends to negative times and to the case of a non-constant a by going back 
to definition (ITTI) of v{t, r). Introduce initial and final values of a{t) as aij, and a = (oi + af)/2. When 
t < 0, a(t) = tti and the previous method straightforwardly provides 

Fit) = - Im Livi + iO+) if t < 0. (73) 
a; 

For t > 0, the jump of a(r) at r = creates a difficulty, in view of the derivative 

du v{t,T)-v{T) ia(r) !T(t, r) - i>(r) i 

d7^'' ^) = t-r + W^) ^ t-r + 2t"^^)- ^^^^ 

Indeed, the discontinuity implies that d(r) oc 5(1"), which we have used in simplifying the last term of (l74l) . 
Thus, a separate treatment of the vicinity of r = in the integral /l^ dr that defines F{t) in Eq. (|40l) is 
required. We write its contribution as (t > 0) 



/e QT 1 
m{vit,T))-diT) 
t — T 2t 

2 /""f / i \ i 2 /""f 

~ - Re / m I t>t -I (of + a) ) — da = - Re / dMm(M), (75) 

t v/oi V 2t J It t Jui 

where Uj = t^f + ia/t and Uf = Vf + iaf/t. Therefore, by m = dp/dv. 



Fr=o{t) = '^Re 

V 



P ( + ly ) - p U'f + 1^ 



(t > 0), (76) 



which is a contribution relative to the momentum transfer induced by core-width variation. In going to the 
second line of Eq. ( 1751 ) an assumption has been used that a(r) varies continuously between values Oi and 
within a time interval shrunk to zero. This amounts to making a constitutive assumption on a{t), indicated 
by the sign that stands for a weak equality in the following sense. Consider some function g{y) of 
v{t) = Vi9{—t) + Vf6{t). Then g{v{t))5{t) has no meaning within Schwartz's theory of distributions, where 
6^(0) is undefined. However, a choice 



g{v{tmt) 





■ 1 








- / dvg{v) 




.Vf — V 





6{t), (77) 



provides the "intuitive" resuk that ^/(t;(t)) = f'{v{t))lvj5{t) = 5{t) J^;.^ dv f'{v) = lf{v)j5{t). In 
general, such weak definitions of the product g5 are not unique (Colombeau, 1985) and the adopted one 
must be warranted by the physical context (Colombeau, 1989; 1990). Here, in a pragmatic way, Eq. (1761 ) 
and its consequence, Eq. (|79l) , will be supported by numerical results in the next Section. 

The remaining part of the integral J'i^ dr is a principal value from which r = is excluded. We split it 
into parts dr and /q+ dr, each one being integrated as in Section[6l but with a(r) = Oi and a(r) = Of, 
respectively. A simple calculation led as in Section [6] yields the following contribution (t > 0)): 
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« 2v{t) 
= wo- 



a(t)cs 



/* dr dv 
m(u(t,r)) — (t,r) 
-oo t — T (XT 



(78) 



Im 



Of 



L h;f + i 



. Of 



where Fr=o{t) shows up without having to mukiply distributions. 
Upon gathering Eqs. ( 1731 ). ( 1761 ) and ( 1781) . F^=o cancels out so that: 



F{t) = - Im L{vi + iO+) if t < 0, 



Im 



Of 



L h;f + i 



. Of 



L(T;f + if) -L(t;i + iO+)" 



bl + if 



if t>0. 



(79) 



This expression reduces to (I63al) for Of = Oi = a and t > 0. The expansion as a^f/t — )■ is identical to Eq. 
(l64l) with a replaced by af. 

Since the initial and final widths Oi ^, as computed from v^f by Eq. (I72al) . are proportional to the interplane 
distance c?, the self-force scales as 



F{t) = wod-^F{cst/d). (80) 

The master curve F, which depends on Vi/cs and Vf/cs as parameters, is such that F{x) ~ F{0^) for 
X <^ 1, and ~ F{oc) + ix~^ for x ^ 1. The latter regime is described by Eq. (|64l) . from which the 
analytical expressions of F(oo) and the constant i are read. One finds for both screw and edge that 

F(0+) = —^lm L{v, + iO+) + 2- f ^ - (81) 
Wo a Cs Vflf a / 

Figs. [3ta) and (b) represent the function F{x) as a function of x > in log/log scale, for screws and 
edges with various terminal velocities. Material parameters a = 0.1 and cl = 2cs are used. The initial 
velocity is Vi = and the terminal velocities are Vf = v, with v/cs = 0.5, 0.7, 0.9, 1.01 and 1.2 for the 
screw (Fig.[3]^a)), and the same velocities with v/cs = 2.1 added, for the edge (Fig.|3]^b)). Expression (I72al) 
of the core width is used. The asymptotic regime given by Eq. (|64l ) is represented for some of the curves 
as dashed lines. The self-force is seen to strongly depend on Vf. A similar conclusion would be drawn upon 
varying vi, which may lead to negative values of F{x) (not shown). 

Fig. ^h) makes conspicuous a wiggled structure for intermediate values of x in the edge case. The 
0(t~^) term in the asymptotic formula (l64l) is seen to be relevant only in a range of high x values in non- 
subsonic cases. The wild variations of the asymptotic regimes are elucidated in Figs.[3];c), (d) where F{0~^) 
and F(oo) (which is independent of Vi) are displayed as functions of the terminal velocity. 

Fig. He) illustrates the dramatic increase of the self-force at v = Cs (screw) and v = cr (edge), as 
the core width shrinks (to a non-zero value due to a). At f = cs, the edge core width blows up and 
the self-force vanishes (see previous section). The quantity F(oo) in Fig. |3]^d) being proportional to the 
radiative dissipation rate, it vanishes at the radiation-free edge velocity Vc = V^Cl (even for a > since 
the phenomenological drag contribution Fd is not part of F). 
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1.5 2.0 

v/cs 

Fig. 3. (a) and (b): master curve F{x) vs. x for (a) screw and (b) edge dislocations subjected to a velocity jump from 
rest to various terminal velocities V{ (see text). Solid: Eq. ( f79l ). The line thickness decreases as Vf increases. Dashed: 
asymptotic expression (l64l ). (c) and (d): asymptotic values of -F(O) for zero initial velocity, and F{oo) (independent 
of the initial velocity) as functions of the terminal velocity v = Vf. A value cl = 2cs is used. 

8.2 Smooth transition in velocity and core width 



The \/d scaling of F(t) reported above makes F{t) blow up at the origin in the Volterra limit d — )■ 0. 
On the other hand and for \v\ < cs, F{t) goes to a finite limit for times such that c^t/d ^ 1, owing to the 
beneficial action of the counter-term examined in Section |6l This blowing-up as d~^ for c^t/d <^ 1 finds its 
origin in the infinite acceleration that takes place at t = 0+ in the velocity jumps that were considered up to 
now. 

To make some progress, behavior under finite accelerations is studied by spreading over a finite time tq 
the velocity transition between Vi and Vf, taking for instance 



= Vit + {Vf -Vi) t- To ( 1 



-t/ro 



(82) 



where accelerated motion begins at position {(0) = 0. The acceleration is now finite, of order ,^(0) 
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{vi — Vi)/TQ. With v{t) = ^{t), we use a{t) = a{v{t)) as given by Eq. (I72al) . The motion reduces to that 
considered in Section [8] as tq — )■ 0. The integral in F{t) over the range r < is done explicitly as before. 
One obtains 



F{t) = -lmL{vi + iO+) for 



t < 



(83a) 




L{v*it))-L{v, + iO+) 

V*(t) — Vi 



for 



t > 



(83b) 



v*{t) 



dcf 



+i[a(t) + ai]/2 



(83c) 



where the remaining integral must be evaluated numerically. These equations hold for any ^{t) at positive 
times, independently of Eq. (l82l) . 

In the following, only velocities \v\ < cs are considered for simplicity, so that F(t) = for t < 0. 

A straightforward series expansion shows that the self-force grows up first linearly with time, as 



This result holds indifferently for screws and edges, and comes out of the explicit contribution in Eq. (|83b|) . 
while the integral itself is an O(t^). 

The dynamical regimes that follow this linear growth step are discussed on the basis of numerical eval- 
uations of F(t) (see AppendixO by means of Eqs. (l82l) and (l83bbc). Figure |4] displays F(t) vs. t > for 
screw and edge dislocations, and two different initial velocities: Vi = 10~'^cs in (a) and (b); and Vi = 0.7cs 
in (c) and (d). In all plots, a = 0.1 (arbitrary value), cl = 2cs, and the final velocity is V{ = 0.9cs, which lies 
slightly below the Rayleigh velocity. Four sets of seven curves are represented in log-log scale in each sub- 
figure. Plots in each set are drawn for exponentially-decreasing values of the interplane distance d = 10~^, 
A; = 0, . . . , 6. Each set uses a different value of the characteristic acceleration time tq: the grey set (dashed) 
uses To = (infinite acceleration) and represents the jump formulas (|79l) ; on the other hand, the orange 
(resp., green, blackiZ] set illustrates the effect of a finite acceleration, with tq = 10 ^ (resp., 10 ^, 10°). A 
marker tags the highest value attainted by F{t). The "edge" plots may possess up to two local maxima due 
to the two intervening sound velocities. Also, the dependence in the initial velocity Vi is less conspicuous 
for the screw than for the edge. In the latter case, F(t) develops a peak for large velocity jumps (Fig. 1Mb)). 
Also, while the highest values of F{t) seem to saturate with l/rf for the screw, a different behavior takes 
place for the edge, whose maximal self-force increases more strongly. 

In all the solid curves of Fig. HI for which tq ^ 0, the initial growth regime with slope one is that given 
by Eq. (f84l) . At large times, all curves collapse to the 1/t asymptotic regime with slope — 1 described by Eq. 
(|64|) with ACT = 0. This asymptote holds for Eq. (|79l) as well, see Section [8TT1 The intermediate regime 
between these initial and asymptotic regimes depends on the acceleration .^(0), as illustrated in Fig.[5l At 
high accelerations, i.e., for values of tq small enough, the self-force reaches at intermediate times a plateau, 
with an undulation in the edge case, of value -^(0"*^) given by Eq. (fSTI) . This plateau is the short-time regime 

^ This order — grey, orange, green, black — applies from the top sets to the bottom ones. See online paper for color 
figures. 




(84) 
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t (a.u.) t (a.u.) 

Fig. 4. Top figures: self-force F{t) vs. time t for (a) screw and (b) edge dislocations subjected to a smooth acceleration 
over a time tq between velocities Vi = 10~^cs and Vf = 0.9cs. The core width varies according to Eq. ( |72a| ). The three 
sets of seven curves correspond to the following values of tq (color online): 5. 10~'^ (black, solid), 4. 10^^ (green, 
solid) and 3. 10"^ (orange, solid). Within each set (from bottom to top), values of the interplane distance d = 10"'^, 
k = 1, . . . , 7 are used. Grey, dashed: asymptotic form (l64l ) with a = Of. Bottom figures: idem, with Vi and Vf as 
indicated. In all plots cl = 2cs and a = 0.1 (arbitrary value). Units are such that cs = 1 but otherwise arbitrary 
("a.u."). 

of the dashed grey curves with tq = 0, and does not depend on acceleration. Such cases are seen in the 
bottom curves of the orange and green sets in Fig. HI 

The plateau begins at a typical crossover time td obtained by matching Eq. (|84|) and (ISTI) . and ends at a 
crossover time tc2 obtained from matching Eq. (|8T| ) and the asymptotic Eq. (|64l ). see Fig. |3a). It is easily 
shown that, for \vi\ < cs, reduces to tq if a is assumed to be velocity-independent. In the general case, 
tci is proportional to tq, unlike tc2 that depends on vi f (both directly and through a^). 

The denomination "high acceleration" means that tq is small enough so that tci < tc2- In the opposite 
case of low acceleration, where t^i > tc2, the plateau is replaced by a mild increase regime, that may end 
up with a peak in the edge case, as illustrated by Fig.|3b). Then, tci,2 are not relevant any more as crossover 
times (see figure). The black curves in Fig. |4] all belong to the latter case. 
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High-acceleration regime 



Low-acceleration regime 



Eq. (81) 





Eq. (81) 





t, 



t 



t 



Fig. 5. Schematic illustrations of the high- (a) and low-acceleration (b) regimes in log-log scale, as determined by the 
ordering of crossover times td and tc2 defined by the intersections of the straight lines. In (a) is reproduced the plot 
of F{t) vs. t of Fig.|4l;b) with tq = 10""^ and d = 1 (orange); in (b), where td and tc2 have no meaning any longer as 
crossover times, data of Fig.Ub) with tq = 1 and d = 10"^ (black) are used. 

The information collected so far for subsonic motion is summarized as follows: 

• For a jump at t = in both velocity and core width between two steady-state regimes: at times t <^ d/cs, 
F scales as 1/d due to infinite acceleration at jump time. In the limit, the singularity generated at t = 
is not Dirac-like, but merely that oc 1/t which arises from the post-jump relaxational asymptote (|67|) . As 
the figures show, this asymptote constitutes an upper bound of F(t) for vi = 0. Thus, except a? jump 
time, the Volterra limit d — )■ of F is non-singular. 

• For finite accelerations: a "quasi-jump" regime is observed, in which the self-force increases in the 
Volterra limit as for jumps. However, this now takes place within a time range of lower bound oc cs/^(0) 
and of upper bound oc d/cs, and occurs only for accelerations and d values that are large enough to pre- 
serve the order of the bounds. By contrast, for low accelerations or for d small enough, a less-than-d"^ 
increase takes place in the intermediate time regime. Its peculiarity is that peaks show up, escaping the 
jump relaxational asymptote, which is indicative of delayed relaxation. 

9 Asymptotic logarithmic behavior 

Provided that velocities are small, the linearized equations of Section |5] are a good approximation, and 
information about motion at finite times under the action of a constant stress applied at t = 0, or of an 
arbitrarily varying stress applied between t = and some instant t = ti is extracted (Eshelby, 1953) from 
Laplace transforms (l58l) : for both screws and edges, the time-dependent quasimomentum reads p{t) ~ 
rrf^{t)v{t), where the function nf^{t) is mo times the logarithm in Eqs. (I55bl) and (I56bl) with \uj\ replaced 
by This is deduced from the small-s expansion of m(s), whose leading term is identical to that in those 
equations, with \uj\ replaced by s. The small-s behavior thus translates into the well-known logarithmic 
time-dependence of the type rrf^it) oc log(t/t*), where t* is a characteristic time proportional to a (i.e., 
to d). The time t* depends on the dislocation character, but also more generally on the loading conditions 
(Eshelby, 1953). As far as the linear theory is concerned, it is remarked that since and are proportional 
to a, the small-s expansion of m coincides with its small-a expansion. Thus, Eqs. (|49|) and (l53l) are only 
logarithmically singular in the Volterra limit. 

As one allows for arbitrary velocities the linear theory is no more relevant, but extracting the exact 
asymptotic behavior of the Volterra limit out of the full-fledged self-force (Ni and Markenscoff, 2008) is 
no easy task. However, one may infer from the latter reference devoted to screw dislocations (and refer- 



26 



a b 




t/t" t/t* 

Fig. 6. Logarithmic behavior of the self-force in the low-acceleration regime, in log-linear scale. Solid (black): plots 
for csTo/d = 10^, 10^, . . . , 10^ (from left to right). Dashed (green): logarithmic asymptotic behavior. Color online. 

ences therein), and from dimensional considerations, that to leading order asymptotic time-behavior of the 
following type should hold: 

F{t)c^m{v{t)) \og{t/t*)m, (85) 

in which t* is proportional to the core width and depends on velocity. 

In the subsonic range, the low-acceleration regime is the place to look for, since the previous Section 
excludes other possibilities. Using the motion of Eq. (|82|) . Fig. |6] displays in the context of Model I log- 
linear plots of the scaled force F{t) /[m{v{t)) ^{t)] vs. t/t* for CsTq/cI = 10^, 10'^, . . . 10®, which ensures 
small acceleration conditions. The constant parameters are a = 0.1 and cl/cs = 2, and the initial and 
final velocities are as indicated. For better reading, the plots do not extend beyond the crossover time td 
introduced in Section [821 The "overshoots" announce the peaks of Fig. IH Prior to plotting the curves, t* 
was determined by plotting the quantity (7 (t) = t exp{—F{t) /[m{v{t))^{t)]} for csTo/d = 10® against cst/d 
in the range (1, 10®) (not shown). It was found that g(t) ~ ci — C2t, where C2 is a very small coefficient, 
of order 10^® ~ 10"''. While this confirms that F{t) behaves as \og{cst/d), it is not clear to this author 
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whether the observed systematic variation of t* with time is a genuine effect or a numerical artifact due to 
precision loss. In practice, t* was identified to ci, which is sufficient to ensure good overall scaling. 

The plots indicate that t* depend on velocity. Moreover, it is observed that even for very small initial 
velocities, the determined values of c^t* /d, indicated in the plots, markedly differ from those provided by 
the leading-order logarithm of the Fourier (or Laplace) mass response function. The latter is read from Eqs. 
(I55bl) and (l57l) . Using the appropriate values of a = a(0) given after Eqs. (I72ab ). this gives c^t* /d ^ 0.540 
(screw) and 1.012 (edge), respectively. To exclude the possibility of an error in the code, F{t) was computed 
for Vi = 0, = 10"'' both by direct numerical integration over time using the full integral formula (see 
AppendixO, and by numerical Laplace inversion of the linearized equation F''"(s) = s'^rn{s)^(s) by means 
of a contour that runs on both sides of the negative-axis cut, skirting round the origin. Same results of F{t) 
were found to four significant digits for screw and edges, which validates both procedures. It ensues that 
the characteristic time provided by the FT or LT is irrelevant to general applications, and that a theoretical 
expression (or other means of computation) of t* that reproduces the values indicated in the plots should be 
looked for. For screws, the expression provided in Ni and Markenscoff (2008) is a possible candidate, but a 
detailed study of this issue lies outside the scope of the present work. 

10 Conclusion 

Starting from exact dynamical core equations for screw and edge dislocations in an isotropic medium, 
equations of motion (EoM) were deduced by means of an arctangent ansatz with arbitrary time variation 
of the core width. The stationary line energy density function W{v) was shown to uniquely determine the 
self-force, with same structure in both cases, Eq. (|40|) . They contain a special local term the regularizing 
nature of which was made conspicuous. This term naturally arises from the formalism employed (Pellegrini, 
2010). It is distinctive of the present theory, although its concealed presence can be traced in the EoM for 
screw dislocations derived by Eshelby (1953) from an electromagnetic analogy. For |f | < cs in the volterra 
limit, this term was shown to cancel out an otherwise diverging contribution in the self-force. An expression 
for it in terms of the energy, Eq. (l66l) . was proposed. 

The arctan ansatz evidently simplifies matters, allowing one to introduce a complex "velocity", of real 
part an average velocity, and of imaginary part an average core width divided by the averaging time in- 
terval, which is perhaps not too surprising in view of the analytic structure of static PN solutions (Lejcek, 
1976). This device, which consistently handles the square-root branch cut, makes the obtained self-force 
expressions formally hold for all velocity regimes, as has been shown for steady motion in Section |71 

Specific analytical results were obtained: the self-force F(t) associated to uniformly moving screw and 
edge dislocations that undergo an instantaneous velocity change, Eq. (|63a| ): its generalization to instanta- 
neous changes in core width, Eq. (1791 ); the frequency-dependent mass and damping function of the edge 
dislocation, Eqs. (I56al) and (I56dl) : Eq. (|69| ). which provides a simple definition of the important functions 
A{v) and B{v) that enter Weertman's equations (1969); the Rosakis (2001) steady-state equations in terms 
of the Lagrangian. Although some of these features are relevant to velocities \v\ > Cs, applications to this 
regime should be regarded with caution for reasons evoked in Section |7] Finally, a synthetic expression of 
F{t) for arbitrary motion from a state of constant velocity was given in Eq. (|83l ). 

To support these findings, various checks the EoM and of the underlying formalism were performed, 
allowing one to retrieve a number of known results as particular cases: in Section [Z2l results by Marken- 
scoff (1980) and Callias and Markenscoff (1980) concerning the self-stress on the slip plane produced by a 
Volterra dislocations; as Eq. (1671 ), the self-force for subsonic smeared-out screw and edge dislocations jump- 
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ing from rest to constant velocity, in the limit of vanishing core width (Clifton and Markenscoff, 1981); in 
Section 1441 Eshelby's EoM for screws (1953); in Section [51 the frequency-dependent mass of the screw 
dislocation in oscillatory motion at small velocities (Pillon et al., 2007); in Section|7]the steady-state EoM 
of Model I (Rosakis, 2001). All of these stand as indirect verifications of the dynamical core equations 
reviewed in Section |2l 

Finally, a distinction was introduced in the subsonic range between high- and low-acceleration regimes 
on the basis of a quantitative criterion, log t behavior of the self-force being observable for low accelerations 
only, associated to a non-trivial characteristic time about which we gave some information of numerical 
nature. No similar analysis was undertaken in the context of sonic transitions, for lack of a proper dynamical 
core- width model. 

Our generic expressions of F{t), which at least apply to screw and edge dislocations in an isotropic 
medium, might extend to anisotropic media as well (BuUough and Bilby, 1954; Teutonico, 1961), although 
this remains in need of a formal demonstration. Obviously, the underlying structure revealed by our analysis 
(e.g., the occurrence of a complex Lagrangian in connection with dissipation in Mach fronts) calls for deeper 
understanding. 
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A Energies and related matters 

The velocity-dependent kinetic energy, Wk, strain energy, W^, and total energy W = Wy, + W^, of the 
displacement field associated to uniformly moving screw and edge dislocations have been computed by 
Weertman (1961), after early steps were taken by Frank (1949) and Eshelby (1949). See Weertman and 
Weertman (1980) and Lothe (1992) for reviews. 



Let /3s,L = Y 1 — (w/cs,l)^- In notations close to that of Hirth et al. (1998) one has, for the screw dislo- 
cation, 





(A. lb) 



(A.la) 



W{v) 



wo/3s , with W{0) = Wo, 



(A.lc) 



and for the edge dislocation 
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2 

W^{v) = ^ (4/3l + 4/3l ^ + /3| - 5/3s - S/^s"' + Pf) , (A.2a) 

W^s(^) = Y (7)' (12/3l + 4/3l ' - /3s' " 9/3s - T/^s"' + /^s"') , (A.2b) 



W^(^) = ^ (^^ j (16/3l + 8/3l ' - 14/3s - 12/3s-' + 2^^) (A.2c) 
with iy(0) = = 2m;o - (A.2d) 



where u is Poisson's ratio. In statics, the quantity wq = log(-R/ro) (/i^^) /(47r) is the rest energy factor. 
The logarithm log(/2/ro), where R and ro are the outer and inner cut-off radii, does not enter the dynamic 
problem (Beltz. et al., 1968), and the identities to be given are independent of it. Accordingly, we use 
everywhere wq as defined by Eq. (|33] ); whenever the terms "energies" or "Lagrangian" are employed, they 
refer to prelogarithmic terms only. 

In the text, we need the limit of W]^{v) as — )■ 00 with Im > 0, and the order of the next term. Because 
with the principal determination of the square root, ^ —z^ = —iz signlm^; for 2 G C \ M, one has in this 
limit /3s,L ~ — i(^/cs,L)- Owing to the term —/3s for the screw, and to the term l^l/v"^ for the edge, the above 
expressions yield for both 

W^k(^)=i^- + 0(1AO, (A.3) 
2 cs 

The "stationary Lagrangian" is (Stroh, 1962; Beltz et al., 1968; Hirth, Zbib and Lothe, 1998) 

L{y) = W^{v) - W,{v) = 2W^{v) - W{v). (A.4) 



Introducing /3s2 = " ^V(2ci), Eqs. SKTkh) and (IXIbb) provide (Bekz et al., 1968) 

L{v) = —wo/3s (screw), (A. 5a) 

= -Awo (^)' (/3l - /3s~Vs2) (edge), (A.5b) 

sothatL(O) = —Wo (screw) and L(0) = —wo/{l — i') (edge). The edge Lagrangian vanishes at the Ray leigh 
velocity where the kinetic and potential energies coincide (Teutonico, 1961). 

A quasimomentum is introduced as p = dL/dv (e.g., Hirth et al., 1998). The following identity was 
demonstrated by Beltz et al. (1968) from the volume-integral expression of L{v): 

Wk = ^pv. (A.6) 

Thus p{v) is readily deduced from the explicit expressions of Wk, and matches Eqs. (|34l ) and (|4TI) . Intro- 
ducing a mass function as m = dp/dv (other definitions are found in the literature, e.g. Weertman, 1961; 
Sakamoto, 1991; Ni and Markenscoff, 2008), it follows that 

dp d'L d /Wk\ , . „, 

m = — = — — = 2— . (A.7) 



dv dv'^ dv \ V 
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From the above, one also sees that (e.g., Hirth et al., 1998) 



^d. ^. IdW ^^^^ 
m = - — {pv- L) = (A.8) 
V av V av 

where the first equality makes use of the definitions p = dL/dv and m = dp/dv, and the second one 
appeals to Eqs. (IA.4I) and (IA.6I) . The right-hand side complies with Frank's (1949) definition of p, such that 
dp/dt = dW/dx (x being the position), which implies in the steady state that 

dvdW ^ ^ 

p= / — . (A.9) 

Jo V dv 

Gathering Eqs. (IA.7D and (IA.8I ) leads to the differential relationship 

2iL m . if:. <A.,0) 

dv \ V / V dv 

Finally, as in relativistic mechanics of mass points (e.g., Bergmann, 1976, p. 92), the kinetic energy of the 
dislocation is the velocity-dependent part of its energy 

f" dW f" 

K{v) = W{v)-W{0)= / dv^{v)= / dvm{v)v. (AM) 

Jo dv Jo 



B Useful integrals 



Some integrals used in the calculations of Section|4]are listed hereafter, where it is assumed that lra{v) > 
0, with Re{v) arbitrary. Principal determinations of the functions are used as defined in Abramowitz and 
Stegun (1972). 



1 du 



u 



1 



-1 TT ^Jl — V?U 

1 du 1 



1 TT ^/I^{u-Vf (1 -y2)3/2' 



1 du 1 



1 TT U'^ 



U 



u^ 



u — V 



du 1 

cs<|m|<cl U^ U — 



Cl 

V 

Cl 



coth 



Cl 



Cl 



coth V 



— - tanh"^-^ \ - { — - tanh"^ — 



cs 

V 
Cs 



coth" 



Cs- 

(0 < Cs < Cl), 

V 
Cs 



(B.la) 
(B.lb) 

(B.lc) 



(B.ld) 



It is recalled (see above reference) that tanh~^2; = coth~^2; + i(vr/2) sign(Im z). 

The following series (Gradshteyn and Ryzhik, 2007) serves to evaluate the Fourier transforms over time 
in Section 151 We use the notation R-iu{z) = ^\Iy{z) — 'hj.lyz)], and Ky, I,, and L,^ are the modified Bessel 
and Struve functions (e.g., Abramowitz and Stegun, 1972). 
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noo e'^^Mt 



(l+t2)5/2 

(1+^2)3/2 

oo gi'^idt 
? 



\u\ iL2(|a;| 



^2)1/2 



|a;|A'i(|a;|) + 1 - ILi{\u\) 
Ko{\uj\) + isign{u)ILo{\uj\), 



rco . , 1 

/ (v^TTt^ - t) e'-Mt = — - ^-llir^ + -iLi(|a;|). 



(B.2a) 
(B.2b) 
(B.2c) 
(B.2d) 



As functions of v in the upper complex plane, integrals (IB.lal) - (IB.ldl) have been checked numerically with 
the Mathematica software (Wolfram Research, 2007). A similar check was carried out for integrals (IB.2al) - 
(IB.2d|) as functions of a; G M. 



C Numerical method 



The integral in Eq. (I83bl) is evaluated with Mathematica' s. general-purpose adaptative integration routine 
(Wolfram Research, 2007), and a change of variable t = t — exp(— m) with u running from — log t to +oo. 
Beforehand, Eqs. (l83l) were compared numerically with the more general Eq. (l40l) for various test cases 
including negative times, and decelerated motion from supersonic to subsonic velocities. Equation (|40l) was 
integrated with u from — oo to +oo with u = — logt specified as a possible numerically problematic point 
when t > 0. The reduction to a finite integration range was left to Mathematica. Same numerical values 
were found, to seven significant digits at least in all cases examined. Prescriptions such as v + iO"*" in Eqs. 
(l72bb) are implemented with 0+ < 10^^. 
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